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Abstract — The purpose of this paper is to improve the 
robustness of traditional image watermarking based on singular 
value decomposition (SVD) by using optimization-based 
quantization on multiple singular values in the wavelet domain. 
In this work, we divide the middle-frequency parts of 
discrete-time wavelet transform (DWT) into several square 
blocks and then use multiple singular value quantizations to 
embed a watermark bit. To minimize the difference between 
original and watermarked singular values, an optimized-quality 
formula is proposed. First, the peak signal-to-noise ratio (PSNR) 
is defined as a performance index in a matrix form. Then, an 
optimized-quality functional that relates the performance index 
to the quantization technique is obtained. Finally, the Lagrange 
Principle is utilized to obtain the optimized-quality formula and 
then the formula is applied to watermarking. Experimental 
results show that the watermarked image can keep a high PSNR 
and achieve better bit-error rate (BER) even when the number 
of coefficients for embedding a watermark bit increases. 

Index Terms —SVD, DWT, PSNR, BER 


I. INTRODUCTION 

With the rapid development of activity on the internet, much 
digital information is widely spread. Digital watermarking 
was developed to hide digital information and protect the 
copyright of multimedia signals, like audio, images, etc. Due 
to the fact that discrete-time wavelet transform (DWT) 
provides a useful platform, numerous DWT-based algorithms 
for digital watermarking have been proposed in recent years. 
Watermarking in the spatial domain [1-11] is usually more 
vulnerable than watermarking in the frequency domain 
[12-29] with the same embedding capacity due to the fact that 
spatial-domain methods are generally fragile to 
image-processing operations and other attacks [23-25]. The 
spatial-domain singular value decomposition (SVD) for 
image watermarking was first introduced by Liu et al. [8]. In 
this paper, the authors used a spread-spectrum technique to 
embed a watermark by modifying the singular values of the 
host image in the spatial domain. Some authors embedded 
watermark to U and V components to increase embedding 
capacity [9, 10] while Ghazy et al. [11] presented a blockby- 
block SVD-based image-watermarking scheme to increase 
embedding capacity. However, the robustness of SVD-based 
image watermarking in the spatial domain is low. In recent 
years, many image-watermarking techniques combine DWT 
and SVD to achieve better transparency and robustness [17, 
18, 24, 25]. Bao et al. [17] proposed a novel, yet simple, 
image-adaptive watermarking scheme for image 
authentication by applying a simple quantization index- 
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modulation process on each single singular value of the 
blocks in the wavelet domain. Their watermarking scheme is 
blind and is robust against JPEG compression but extremely 
sensitive to malicious manipulation such as filtering and 
random noising. Ganic et al. [18] applied SVD to all details, 
approximating part of the DWT and watermark image to 
increase embedding capacity. Gaurav and Balasubramanian 
[24] embedded a watermark into the reference image by 
modifying the singular value of the reference image using the 
singular values of the watermark. The robustness is slightly 
enhanced. However, the computation is significantly 
increased. Lai and Tsai [25] reduced the computation in [24] 
by directly embedding the watermark into the singular values 
in the wavelet domain. In this work, we first divide the DWT 
middle frequency parts LH3 and HL3 into several square 
blocks to have high embedding capacity. Unlike the 
traditional spread-spectrum technique on single singular 
values [24, 25], we use multiple singular value quantizations 
to embed a watermark bit. It does not only keep a high 
embedding capacity but also achieves strong robustness 
against median filtering. On the other hand, an optimized 
quality formula is proposed by minimizing the difference 
between original and watermarked singular values. First, the 
peak signal-to-noise ratio (PSNR) is defined as a performance 
index in matrix form. Then, an optimized quality functional 
that relates the performance index to the quantization 
technique is obtained. Finally, the Lagrange Principle is 
utilized to obtain the optimized quality formula; then, the 
formula is applied to watermarking. Experimental results 
show that the watermarked image can keep a high PSNR and 
achieve a better bit error rate (BER) even when the number of 
coefficients for embedding a watermark bit increases. In 
particular, the robustness against median filtering is 
significantly improved. 

This paper is organized as follows. In Section II, we review 
some mathematical preliminaries. Section III introduces the 
proposed watermark embedding and extraction. In Section 
IV, we rewrite PSNR as a performance index. An 
optimized-quality equation that relates the performance index 
to the quantization constraint is proposed, and the Lagrange 
Principle is used to solve the optimized-quality problem. The 
solution is utilized to embed the watermark, and we discover a 
very good result; the watermark is extracted without the 
original image. In Section V, we present some experiments to 
test the performance of the proposed scheme. Finally, 
conclusions are drawn in Section VI. 


II. PRELIMINARIES 

In this section, some related steps for the proposed image 
watermarking scheme are reviewed. 
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DISCRETE-TIME WAVELET TRANSFORM (DWT) 

The wavelet transform is obtained by a single prototype 
function which is regulated with a scaling parameter and shift 
parameter [28-31]. The discrete normalized scaling and 
wavelet basis function are defined as follows: 

p,-,/( t ) = 2 'M 2/( - r ) (1) 

( 2 ) 


where j and x are the dilation and translation parameters; 
from this, one can require that the sequence 


{0}c---cV' 1 c\/ o c\/_ 1 c---cL 2 (R) 

(2) 

Forms a mutiresolution analysis of L2(M) and that the 
subspaces • Wl, WO, W_ i, • • • are the orthogonal 
differences of the above sequence; that is, Wj is the 
orthogonal complement of Vj inside the subspace Vj - 1. 
Then, the orthogonality relations follow from the existence of 
sequences h = {hx}x E 1L and g = {gx}x E H that satisfy the 
following identities: 

h T = Po,o.?-i,r and p(f) = 'J2j^h T <p{2t-T) 


T€Z 


(4) 


g T = 00,0,p-i, r and 

(5) 

where h = {hx}x E Z and g = {gx}x E Z are, respectively, the 
sequence of low-pass and high-pass filters. In this paper, we 
use a Haar scaling function and wavelet to transform the host 
image into the orthogonal DWT domain by three-level 
decomposition. A method to implement DWT is a filter bank 
that provides perfect reconstruction. DWT has local analysis 
of frequency in the space and time domains, and it obtains 
image multi-scale details step by step. If the scale becomes 
smaller, every part gets more accurate and ultimately all 
image details can be focalized accurately. If DWT is applied 
to an image, it will produce high-frequency parts, 
middle-frequency parts, and a lowest-frequency part. Figure 1 
shows the procedure of applying one-level DWT to an image. 
In order to guarantee both image quality and robustness, this 
study embeds the watermark into the middle frequency parts 
LH3 and HL3 in DWT level-three. 


SINGULAR VALUE DECOMPOSITION (SVD) 

The singular value decomposition of a matrix A with size m x 
n is given by 


A = UDV r 


( 6 ) 

Where, U and V are orthogonal matrices, and D = diag(Ai) is a 
diagonal matrix of singular values Xi, i = 1,2, •••, which are 
arranged in decreasing order. The columns of U are the left 
singular vectors, and the columns of V are the right singular 
vectors of image A. 



Figure 1 2D DWT 

III. OPTIMIZATION SOLVER 

To find the extreme of the matrix function, some optimization 
methods are summarized in [29-31]. The operations of the 
matrix function are first shown as follows: 


Theorem 1: If W is a k x k constant matrix, and X A is a k x 1 
column vector with k unknown variables, then 


( 7 ) 

Theorem 2: If X is a k x 1 constant vector and A X is a k x 1 
column vector with k unknown variables, then 





Hast Image 



Watermarked image 


Figure 2 Watermark Embedded Process 
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(13) 


Now we consider the problem of minimizing (or maximizing) 
the matrix function f (X A )subject to a constraint 
b XVa 0. This problem can be described as follows: 


minimize 



9(a) 


subject to gf X] — 0 

X ' 9(b) 

Theorem 3: Suppose that g is a continuously differentiable 
function of X on a subset of the domain of a function f. Then 
if Xo minimizes (or maximizes) f (X)subject to the Constraint 
(X) ; Vf (Xo) and ; Vg (So) are parallel. 


That is, if b_X0_^0, then there exists a scalar such that 

(10) 


Watermarked image 



Watermark B 


V_ 

Figure 3 Watermark extraction process 

"(*-*)='(*)+*(*) 

Then the original problem (9) becomes a function H (X£) 
which has no constraint. The necessary conditions for 
existence of the extreme of function H (X^) are: 



( 12 ) 


IV. Proposed Optimization-based DWT-SVD 
WATERMARKING SCHEME 

The proposed watermarking scheme is introduced in this 
section. The watermark is extracted without the original 
image. 

WATERMARK EMBEDDING 

The proposed embedding process as shown in Fig. 2 is 
summarized as follows. 

(1) Use three-level DWT to decompose the original image A 
into four subbands (i.e., LL3, LH3, HL3, and HH3). 

(2) Divide LH3 and HL3 into non-overlapping blocks Ak, k = 
1,2, *-,N. 

(3) Apply SVD to each block, i.e., 

A k _ ukjyk/■yk\ T k = L 2 . _ N 

(14) 

where k represents the number of blocks in LH3 and 
HL3. 

(4) Watermark B = {|3j} randomly generated using a 
binary PN sequence is embedded by modifying 



Figure 4 Optimal Embedding 

Singular values X X i = 1, • ••, r of the matrix A K as follows: 
Let 
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WATERMARK EXTRACTION 


(15) 


To detect the watermark, DWT is first performed and then the 
optimized singular values of SVD in each block are grouped. 
The embedded binary bits are extracted by using the 
following rule. 






Figure 3 shows the detailed process of the proposed 
watermark extraction. 


V. OPTIMIZATION OF PSNR ON SINGULAR VALUES 


Generally, the quality of a watermarked image is evaluated by 
the peak signal-to-noise ratio (PSNR). Since a tradeoff exists 
between image quality measured by PSNR and robustness 
measured by BER, a scalar parameter £, is applied to connect 
the PSNR and the quantization equation to optimize the 
tradeoff in this section. The details are in the following: 



Figure 5 The relation b/w PSNR and number of coefficients 


TABLE 1 


The PSNR. If I(i, j) and I(i, j) are the values of the original and 
the corresponding modified pixel in the original image I and 
watermarked image I, then is defined as 


PSNR — Ifljfqp 


10 


25 S 2 phpj 


16 


where m and n represent the height and width of the host 
image. Based on the watermark embedding of singular values, 
PSNR is expressed as 


PSNR - -1 tilog 


m 


255 2 j 


(17) 


For the optimization of the watermarked image quality, 


Eq. 17 is rewritten as a performance index: 

||V“-V || £ 


/(V) = J 


255 2 wijj 


Or, 


f/l \ “Aj/: J j 

AM 


(18) 


(19) 


-A^y j —Ajy 

25 5 2 mti 


mirumizc- 


£.u b j c c t to WXpt — y 1 


20(a) 

20(b) 


To embed the watermark B, we need to solve the optimization 
problem (20). By Theorem 3, we set the Lagrange multiplier X 
to combine (20a) and (20b) into a matrix function: 


H 




( 21 ) 


which has no constraint. Since 2552 mn is a constant, we 
redefine H_VDk ; ^_as follows: 

H ^A|y , ^ J — ^A^y-A^J ^A^y-A^yJ 
+ 25S 2 wm£(wA Dt -y 1 ) 

2 E - 2 ( V-A I y ) + 25 S 2 ™*fW T - 0 
3.A jy ■ ' 

"jLT 


( 22 ) 

23(a) 
23(b) 

(24) 

Since W A XDk !4 yl from (23b) and 255 2 mn is a scalar, we 
rewrite (24) as 

2SS 2 mii 


Multiplying (23a) by W, we observe that 
J^WA^-WA^) + 25S 2 »!«£WW T - 0 


(^-WA^) 


-(TO 1 - 0 


(25) 

Some operations yield the optimal solution for parameter 
as 

(26) 

Replacing Eq. 26 with Eq. 23a yields the optimal embedded 
singular values 


= Atf-W’ (WWT (WV-J.) 

By using y2 instead of yl yields the optimal embedded 
singular values 

V* = A rf -W T (WW T ) "‘(WA^-Vj) 


(27) 


(28) 
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VI. RESULS 

This section presents experimental results that indicate the 
performance of the proposed image-watermarking scheme. 
Forty host images including the four images, Lena, Jet, 
Peppers, and Cameraman, each a size of 512 x 512, are 
decomposed into three levels by applying DWT; then, the 
watermark is embedded into the LH3 and HL3 coefficients. 
Figure 5 shows that the watermarked image can keep a high 
and stable PSNR (almost 54.5 dB) even when the number of 
coefficients for embedding a watermark bit increases. This 
feature indicates the proposed optimization embedding 
formula using Lagrange principle. In order to compare with 
the SVD-based method [25], PSNR is fixed to be 55 dB. 
Table 1 shows the comparison of the embedding capacity 
under fixed PSNR = 55. Figure 6 shows the original images, 
and Figs. 7 and 8 show the watermarked images obtained with 
different parameters. 



Figure 6 Original images 



lal Lm tf=55i cl* i Wi tpMi 





Figure 7 Watermarked image for k=4 



Figure 8 Watermarked image for k=8 


To evaluate the robustness of the proposed method, the 40 
images including the four images, Lena, Jet, Peppers, 
Cameraman, are tested. After the embedding process, four 
attacks are adopted to test the robustness of the embedded 
watermark in cases k = 4 and k = 8. The robustness is 
measured by BER defined by 

BEE - 100% 

Where, 

B error and B total denote the number of error bits and the number 
of total bits, respectively. The method proposed herein is 
compared with the SVD based method using the 
spread-spectrum technique in the DWT domain [25]. 

VII. Conclusion 

This study improved the robustness of traditional SVD based 
image watermarking by using optimization-based 
quantization on multiple singular values in the wavelet 
domain. Experimental results show that the watermarked 
image can keep a high PSNR and achieve a better BER even 
when the number of coefficients for embedding a watermark 
bit increases. In particular, the robustness against JPEG 
compression, Gaussian noise, and median filtering is 
significantly improved. The future work is the consideration 
of improving robustness against rotation. 
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